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FIG.29 



rea J*4 ep (5> . x C51 , y (5) , z (5) , xs (5) , ys <£) , zs (5) 
rea 1*4 al(20CKbl (200) . ol C20t) 

iij=1 234556 
f0=0. 05 

dlh=15.0/f0 [ cm 
na=16 
write (*. 10) 

forwatC Enter the location of x. y,z (cm) : ' $) (2) 
read{*. *, end=90> xp, yp, zp 



caM marrayCxp. yp,zp,na, cl) (3) 
do i=2,na+1 

verr=ran<i i j) 

at <i-1)=cl 0^Cverr-0. 5)*0> 001)-cl (1) 

end do 

write<*.*) AL(cm)\ (al (i>, i=1,na) 
wr ite(*,*) 

call ncycla(na, dlh, at) (5) 

do 1=1,5 

ep(t)=1.0e20 
and do 

do tx=-30, 30 

xp=f faat(ix)*J0. 0 
do iy=— 30, 30 

yp=f loat (iy)*10. 0 
do fz=-30.30 

zp=fioat(iz)*10. 0 



! noisa 0. 



(4) 



] 



(6) 



(7) 



ca 1 1 aiarr ay (xp, yp f zp. na, c I ) 
do i=2. na+1 

bl<M)=cl(l)-cl<i)-al<M) 
end do 

oa!l mcycle(na,dln, bl) (9) 
er=0, 0 
do i=l. na 

er=er+blC0**2 (10) 
end do 
do i=1,5 

if (er . It. ap(i)) than 
if (f .na. 5) then 
do j=5. 

ep(j)=«p(j r i) 

y<i)=y(j-l) 

2<j)=Z(j-1> 

end do 
end If 

x(i)~xp 

y 0)=yi> 
z(D=zp 
go to 30 
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end if 
end do 
cont i nue 
end do 
end do 
end do 

do 5 

xs(i)*x(i) 
ys(i)=7(i) 

end do 

write(*,*) ' RMS error (cm) x y 

do S=U 5 

write<*,*) sqrt(ep(t)/fIoat(na)), x(i), y (i).z(i) 



do m=1,5 

xO=X8 (m) 
yO=ys(m) 
20=25 (r&) 
do ix=-15. 15 

xp=f !oat Cix)+xG 
do iy=-15, 15 

yp=f I oat (iy) +y0 
do iz=-15. 15 

zp=f (oat (iz)+z0 



call marray (xp,yp.zp, na. of) 
do i=2, na+1 

bl(i-l)=c! (i)-cf (l)-af 
end do 

call mcycle(na,dlh,bl) 
er=0. 0 
do na 

er=er+b! (i)**2 
end do 

do t=1. 5 

if (er . It. ep(i>> then 
if U .ne, 5) then 
do j=5. 

ep(j)=sp(j-1) 



30 



and do 




x (j)=x 



end if 



ep(i) =er 




end if 



go to 35 



end do 
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FIG.31 

cant i nue 
end da 
end do 
end do 

end do 
wr i te {*, *) 

wr ite(* r *) sqrt (ep<1) /float (na» . x (1), y (1) , z (1) (14) 

wr i te (*> *} 
go to 20 

stop 
end 

subroutine marray (xp, yp, zp. na, ci) 
resi*4 ci (1) 

of (1)=sqrt(xp^p+yp*yp+(zp+50. Q)**2) 
do i-2. na+1 

f xx=i/3 

i yy= i - i xx*3 

xm=f I oat { i xx- 1 ) *50. 0-1 0. 0 
ym=f I oat (iyy-1)*50, 0+10. 0 
ot ( i ) =sqrt «xp-xm) **2+ (yp-ym) **2+zp*zp) 
end do 

return 
end 

subroutine mcycle (na, dlh\ al) 
real*4 ai O) 

do i=1,na 
continue 

if (al (i) .gt. dfh) then 
al (i)=a! (i)-dih 
if (al(i) . !e. dih} go to 46 
go to 40 

end if 

continue 

if (al(i) , It. -dlh) then 
al(i)=at(i)+dlh 
if (ai{i) ,ge. -dfh) go to 46 
go to 45 

end if 

continue 
end do 

return 
end 
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FIG.32 



Enter the location of x, y, z (cm) : 152,-203,56 
^ C S7L< 67 - 67562 -38.21133 -1.487458 39 09471 

097327?! 'Vn^ll 11- l^i " 9 " 74537 

-51 7 46763 :3°2693l6 ^ ^ ~ 5 570741 ' 98 28325 



RMS error 
0. 6834297 
0. 8562734 
1. 116775 
1. 163736 
1. 216863 



(cm) 



150. 0000 
150. 0000 
1 50. 0000 
1 60. 0000 
160. 0000 



8.4395386E-02 152. 0000 

Enter the location of x, y, 
ZlL (cm) 5. 506481 

55.74849 14.06553 
-2. 332703 89. 04320 

90.37129 69.39222 



RMS error (cm) 
1 . 445567 
1 . 754374 
1. 951296 
2. 345274 
2. 709345 



-20O. 0000 
-190. 0000 
-20O. 0000 
-230. 0000 
-220. 0000 

-203. 0000 



60. 00000 
50. 00000 
50. 00000 
70. 00000 
60. 00000 

56. 00000 



z (cm) : 22, 123, -89 

57.46710 16.50204 -17 27929 

-20.41722 66.89948 28 19106 

55. 22502 29. 83902 119. 4193 



20. 00000 
20, 0O000 
20. 00000 
20. 00000 
20. 00000 



6. 2024966E-O2 22. 00000 

Enter the location of x, y, 
ZJL(cm) -23.45399 

21. 66080 -27. 96993 

-74. 96463 36. 05470 

18.86572 -15. 62937 



y 

130. 0000 
130. 0000 
120. 0000 
120. 0000 
140.0000 

123. 0000 



-90. 00000 
-100. 0000 
-80. 00000 
-90. 00000 
-100.0000 

-89. 00000 



z (cm) : 60, 161, 5 

32. 54938 -1 3. 85323 -57. 41031 

-77. 36571 22. 85288 -26. 38201 

-9. 367880 -51. 50449 59 00156 



RMS error 
1. 358104 
1 . 400364 
1. 561480 
1. 779230 
1. 850774 



(cm) 



60. 00000 
60. 00000 
60. 00000 
60. 00000 
60. 000O0 



160. 0000 
160. 0000 
170.0000 
1 70. 0000 
150. 0000 



4. 4650473E-02 60. 00000 161. 0000 

Enter the location of x, y, z (cm) 



10. OOOOO 
0. 000O00OE+0O 
0. 0O0O00OE+0O 
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5. 000000 



